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FLUID LIMITS FOR PROCESSOR SHARING QUEUES 
WITH IMPATIENCE 

H. CHRISTIAN GROMOLL*, PHILIPPE ROBERT, AND BERT ZWARTt 



Abstract. We investigate a processor sharing queue with renewal arrivals and 
generally distributed service times. Impatient jobs may abandon the queue, 
or renege, before completing service. The random time representing a job's 
patience has a general distribution and may be dependent on its initial service 
time requirement. A scaling procedure that gives rise to a fluid model with 
nontrivial yet tractable steady state behavior is presented. This fluid model 
model captures many essential features of the underlying stochastic model, 
and it is used to analyze the impact of impatience in processor sharing queues. 
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1. Introduction 

Processor-Sharing Policy and Impatience. Processor Sharing (PS) policies 
were originally proposed as models of time sharing in computer operating systems. 
Recently, generalizations of this discipline have been used to describe data transfers 
in congested routes through the Internet, see Roberts and Massoulie J2H] and Kelly 
and Williams JHj and the references therein. This has created considerable renewed 
interest in the analysis of PS policies. 

This paper studies the behavior of a GI/GI/1 queue serving impatient jobs 
according to the PS policy: if there are N jobs in the queue, each job receives 
simultaneous service at rate 1/N. An impatient job has a random initial lead time 
in addition to its service time. Such a job has a deadline equal to its arrival time 
plus its initial lead time; if the job has not completed service when the deadline 
expires, it abandons the queue (or reneges) and therefore does not complete service. 
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2 H.C. GROMOLL, PH. ROBERT, AND B. ZWART 

For example, the timeout of a TCP flow through the Internet can be thought of as 
the expiration of a random deadline and subsequent reneging of the flow. 

The impact of impatience on PS queues is larger than for First In First Out 
(FIFO) queues. A typical job that abandons a FIFO queue will do so while waiting 
to begin service. In contrast, a job that abandons a PS queue will have already 
received partial service. Since this partial service is wasted, impatience may create 
a significant overhead for a PS server. 

There is a large literature on queueing models with impatience under the FIFO 
discipline. An early paper by Barrer pQ considers an example arising in a military 
application. Stanford |31| is a survey of the literature in this domain (see also 
Stanford [23] and Boots and Tijms [S])- This body of work focuses primarily on 
exact performance analysis. Ward and Glynn {^Sj have recently obtained a diffusion 
approximation for single channel queues. There are also various studies of multi- 
server queues with abandonments, motivated by call center applications; see the 
survey by Gans et al. |10j and references therein. 

There is some related literature treating other policies, but in the context of 
soft deadlines. Jobs with soft deadlines are not impatient; they remain in the 
system until completing service, even if their deadlines have expired. In particular, 
these queues are work conserving. Results for such models describe the extent to 
which overdue jobs are produced by the underlying service discipline, without the 
effect of abandonments. Doytchinov et al. [!5], Kruk et al. [2U1[^, and Yeung and 
Lehoczky |34j investigate the heavy traffic behavior of various systems using the 
Earliest Deadline First and FIFO policies. Gromoll and Kruk J21 describes the 
heavy traffic behavior of a PS queue incorporating a fairly general structure of soft 
deadlines. 

For PS queues with impatience however, only a few results are known. Coffman 
et al. [xj cover the special case of exponential service times and lead times, where 
the lead time and service time are independent. Guillemin et al. |14j consider 
heavy tailed service times, and obtain some results on the reneging behavior of 
large jobs by analyzing the tail behavior of the sojourn time distribution. Using 
some approximations, Bonald and Roberts [J analyze the steady state of a system 
with general service times and some dependence between service times and lead 
times. 



Results of the Paper. This paper analyzes the PS queue with impatience by 
using fluid limits. The dynamics of the system are represented as a measure valued 
process: the system state at time t > is represented by a random point measure 
Z(t) on (0, oo] x (0, oo], such that Z{t) has a point mass at (b, d) € (0, oo] x (0, oo] 
if and only if there is a job in the system at time t with residual service time b and 
residual lead time d. See Jean-Marie and Robert JH] and Doytchinov et al. |§] for 
an analogue representation of residual service times in single server queues. This 
setup enables a fairly general analysis. The case of a general joint distribution of 
service times and initial lead times, with possible dependence of the two random 
variables is included in our setting. 

Under mild assumptions, it is shown that, with a convenient scaling, a family of 
measure valued processes associated with (Z(t)) is tight and converges in distribu- 
tion to some (£(£)). For t > 0, £(£) is a nonnegative measure on (0, oo] x (0, oo], it 
is the limit in distribution of the sequence of random points describing the queue. 
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This fluid limit is characterized as the solution of a functional Equation lj2.8|) which 
can be viewed as a time changed functional differential equation. 

The overloaded case p > 1, which forms our main focus, presents a nontrivial 
and quite interesting steady state behavior. The total fluid mass in the system 
at equilibrium (the fluid analogue of the total number of jobs) is shown to be the 
solution Zqo of a simple fixed point equation 13.2(1 . Moreover, the fluid steady state, 
i.e. the limit of Q(t) as t goes to infinity, is a distribution on (0, oo] x (0, oo] which 
has a simple expression (I2.II|) in terms of z^ . 

These results give also a significant insight on the qualitative properties of PS 
queues with impatience. An interpretation of the fixed point equation l|3.2[) is given 
and used to analyze the total number of jobs in the system and to estimate the 
fraction of jobs that renege. The impact of the variability of the service times and 
of the lead times and other properties of this queue are extensively investigated in 
Gromoll et al. |11| . 

In contrast to the models studied previously in this domain, the service discipline 
considered here is not work conserving. For this reason, analysis of the fluid model 
is more intricate. This is an important difference from earlier work on standard 
PS queues where the fact that the workload process coincides with that of FIFO 
discipline was a crucial ingredient in the proof of the key results. A different 
approach to prove existence, uniqueness, and convergence to steady state of fluid 
model solutions is proposed. It is shown that there exists a maximal fluid model 
solution and by using monotonicity arguments, the properties of the fluid limits 
can be investigated under quite general assumptions. 

Organization. The paper is organized as follows. A detailed description of the 
model and the main results is presented in Section 2. Qualitative properties of the 
fluid model are analyzed in Section 3. Section 4 is devoted to examples. Section 5 
and 6 are concerned with convergence towards the fluid limit. Section 5 establishes 
tightness, and Section 6 characterizes limit points. 

2. Model Description and Results 

This section gives a detailed description of the stochastic processes associated to 
this queue and a summary of the main results. 

2.1. Stochastic model. The stochastic model consists of the following: a proces- 
sor sharing server working at unit rate from an infinite capacity buffer, a collection 
of stochastic primitives E(-), {Bi, D i \'^L l describing respectively the process of ar- 
rivals and the services and the deadlines of the customers, and a random initial 
condition specifying the state of the system at time 0. All random objects are 
defined on a probability space (fi,^",P) with expectation operator E(-). 

The exogenous arrival process (E(t), t > 0) has rate A>0, it is a delayed renewal 
process starting from zero, with ith jump time U. For t > 0, E(t) is the number 
of jobs that arrive to the buffer during (0, i\. For i > 1, Ui is the arrival time of job 
j; jobs already in the buffer at time are called initial jobs. 

For i > 1 , the service time Bi is a strictly positive random variable representing 
the amount of processing time that job i requires from the server. The random 
variable Di is strictly positive and determines the deadline of job i: it represents 
the maximum amount of time that job i will stay in the buffer. Since it arrives at 
time Ui, its deadline is at time Ui + Di. It will abandon the system at this time if 



4 H.C. GROMOLL, PH. ROBERT, AND B. ZWART 

it has not yet completed service. The random variable Di is called the initial lead 
time of job i. 

The model allows either the service time or the initial lead time (but not both) 
to be equal to infinity. Therefore, the random variable {Bi,D{) has values in the 
space R + = [0, oo] x [0, oo]. Here, R + = [0, oo] is the usual compactification of 1Z+ 
with the arithmetic extensions x + oo = oo for all x G M+ , x ■ oo = oo for x > 

2 

and • oo = 0. The collection of Borel subsets of R + is denoted by £S. Throughout 
the paper, it is assumed that all sequences of services and deadlines {Bi,Di}°^ 1 are 

2 

independent and identically distributed (i.i.d.) R + -valued random variables, and 

2 

that their common joint distribution $ on K + satisfies 

#({0} x I) = #(! x {0}) = tf((oo, oo)) = 0. 
Note that the random variables Bi and Di may be dependent. A generic random 

2 

element of R + with distribution ■d will be denoted (B,D). Let p = XE[B] denote 
the traffic intensity of the system. It is assumed throughout that p > 1, that 
is, the server is nominally overloaded. In this way, the classical PS queue, the 
infinite server queue, and mixtures of the two are special cases of this model. (For 
example the GI/GI/oo queue corresponds to the case when service times are equal 
to infinity.) 

Initial condition. The initial condition specifics Z(0), the number of initial jobs 
present in the buffer at time zero, as well as the service times and initial lead 
times of these initial jobs. Assume that Z(0) is a non-negative, integer valued 
random variable. The service times and initial lead times for initial jobs are the 
first Z(0) elements of a sequence (B^,D^) of i.i.d. random variables taking values 

2 

in {(0, oo] x (0, oo]} \ (oo, oo) almost surely. A generic random clement of M + 
distributed as (B®,Dq) will be denoted by (B°,D°). Assume that the expected 
number of initial jobs is finite: E[Z(0)] < oo. 

Time Evolution of the Queue. For each t > 0, let Z(t) denote the number 
of jobs in the buffer (or queue length) at time t, and S(t) denotes the cumulative 
service time per job provided by the server up to time t. Because of the processor 
sharing policy, the quantity S(t) is given by 

(2-1) S(t)= f -L- da, 

Jo Z(s) 

where the integrand is defined to be zero when the queue length equals zero. If a 
job arrived at time s > and is still present in the queue at t > s, at time t it has 
received the cumulative amount of processing time S(t) — S(s). 

Therefore, the residual service time at time t of job i < E(t) (and of initial job 
j < Z(0)) are given by 

Bi(t) = (Bi - (S(t) - S(Ui))) + , and B?(t) = (B? - S(t)) + . 

Define the lead time at time t of job i < E(t) (and of initial job j < Z(0)) by 

(2.2) Di(t) = (Ui + Di-t) + , and £>° (t) = (D® - 1)+ . 

A job's residual service time is the remaining amount of processing time required to 
fulfill its service requirement; its lead time is the remaining time until its deadline. 
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Figure 1. Dynamics of the measure valued process Z(-) 



Job i will depart the system cither when its service requirement is fulfilled or when 
its deadline arrives, it will leave the system at time 

inf{£ > Ui : mm{B t (t), D t (t)} = 0}. 

The state descriptor is a measure valued process that keeps track of the residual 
service times and lead times of all jobs in the buffer. For job i, this information is 

2 

represented as a unit of mass at the point (Bi(t),Di(tj) G R + at all times t > Ui 
such that job % is still in the system. Let 57, 1 denote the Dirac point measure at 

(x,y) s K + if min{x,y} > 0, otherwise 5f ■. is the zero measure. Then the state 
of the system at time t > is represented by the random point measure 

Z(0) E(t) 

Z^ 5 (B«(t),D«{t)) + 2^ 6 (B t (t),D z (t))- 
3=1 t=l 



Z(t) 



Note that the queue length at time t is given by the total mass of the measure Z(t), 

Z(t) = (l,Z(t)), 

where (/, /i) = J-g2 fdfi for a Borel measure /i on R + and a /i-integrable function 

/:!+-► M. 

In this way the dynamics of the system are represented as a distribution of point 

2 

masses on R + moving toward the axes. At time t > 0, points move left at rate 
1/Z(t) and down at rate 1. (A point with one coordinate equal to infinity will 
remain that way while the other coordinate moves.) Point masses vanish when 
hitting one of the axes: a point mass reaching the vertical axis corresponds to a job 
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completing service, while a point mass hitting the horizontal axis represents a job 

abandoning the queue. See Figure ^ 

— 2 
Let Mi denote the space of finite non-negative Borel measures on R + , endowed 

with the topology of weak convergence: ( n — > C m Mi if and only if (/, ( n ) — » (/, £) 

2 2 

for all continuous functions / : R + — > R (recall that R + is compact for the induced 
topology). Let D([0, oo),Mi) denote the space of cadlag paths in Mi, endowed 
with the Skorohod Ji-topology. Then, for t > 0, Z(t) is a random element of Mi 
for each t > 0, and i?(-) is a random element of D([0, oo), Mi). 

It is clear that, given stochastic primitives (E(-), {Bi, Di}^^ and the initial 
condition 2(0), the equation (|2.1(l uniquely determines the processes S(-), Z(-), 
Z(-), and the residual service times and lead times. It is also easily seen that the 
state descriptor Z(-) satisfies the following equation: for each Borel set A E 3&, and 
all t > 0, 

E(t) 

(2.3) Z{t)(A) = Z(0) (A + (S(t),t)) + J2 1 a (Mt), A(*)) , 

8=1 

where A + w — {a + w : a E A} and l\(w) = (1^,^+). Note that the quantity 
Z(Q) (A+ (S(t),t)) corresponds to a shift by the quantity (S(i),t) of the initial 
points: indeed, if (x,y) 6 R^_ and (s,t) 6 R^_, for Aef, 

5(x,y)(A+ (s,t)) = 8( x _ s>y _ t )(A). 

This equation plays a crucial role in determining fluid limits for the model. 

2.2. A Fluid Scaling. A sequence of renormalized stochastic processes (Z (t)) 
associated to the solution of the evolution equation l|2.3|l is introduced. The limits 
of (Z (£)) will give the fluid limits of this queue. 

Let 1Z C [0, oo) be a sequence increasing to infinity. Suppose that for each r S 1Z, 
there is a stochastic model as defined in Section l2~Tl That is, for each r € 1Z, there 
are stochastic primitives (E r (-), {BT, -D[}^i) with associated data A r and f? r , and 
an initial condition Z r (0) which give stochastic processes Z r (-), S r (-), Z r (-), and 
residual service times and lead times {£[(•),£)£(•)} and {B^ r (-), D® r (-)}. Each 
model is defined on a probability space (fT, J? r , P r ) with expectation operator 

£''(•)■ 

A fluid scaling is applied to each model in the sequence. To obtain non-trivial 
scaling limits, initial lead times {D?} will be assumed to be of order r. For each 
r e 1Z, let f? 7 " E Mi be the probability measure defined by 

d r {F x G) = d r (FxrG), 

for all Borelian subsets F,G E S3 with the notation rG = {r ■ g : g E G}. Note 
that if {-BJ",-D[} = {Bi,rDi} for some sequence {i?i,D,}, then ■d r is simply the 
distribution of (Bi, Di). 

For each r E 1Z, the fluid scaled state descriptor is defined, for t > 0, as the 
random measure Z (t) £ M such that 

Z r (t)(F xG) = iz r (rt)(i^ x rG), 
r 

for all Borelian subsets F,G € £%. This definition scales lead times by a factor r _1 
as well. Fluid scaled versions of the remaining processes are defined as follows: for 
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all r £ 11, t > s > 0, and i = 1, . . . , E r {rt), let 

E r (i) = -E r (rt), Z r (t) = -Z r {rt), 

r r 

~S r (t) = S r {rt), Iffat) = 5 r (i)-5 r (s), 

mt) = sr(rt), sT( t ) = I^( r i). 

The following asymptotic assumptions are needed. Let (A, ■&, Co) be fluid model 
data satisfying the assumptions of Section 12.31 Assume that as r — > oo, 

(2.4) ^(O-AO, 

(2.5) tf r -^tf, 

(2.6) Z (0) — ► Co i in distribution, 
in particular, Assumption l|2.4|l implies that A'' — ► A holds. 

2.3. Fluid model. A deterministic fluid model satisfying dynamic equations anal- 
ogous to H2.3fl is introduced. It will be shown later that these equations can be 
obtained as limits of Equation l|2.3[) under an appropriate scaling procedure. 

Let p > 1, and Co G Mj be a measure on K + such that the projections Co(' x K+) 
and Co(K+ x ■) are free of atoms in [0, oo) and z = Co(0, 0) is the total mass of Co- 
Definition 2.1. A measure valued fluid model solution for the data (A,$, Co) * s a 
continuous function £(•) : [0, oo) — » Mi such that 

(i) inf i>a z(t) > for all a > 0, 
(ii) for all C £ C and t > 0, 

(2.7) C(*)(C) - Co(C + (5(0, t),t))+X f 0(C + (5(s, t), t - *)) da, 

Jo 

where S(u,v) = J l/z(s)ds for all v > u > and z(-) is i/ie total mass 
function is z(-) — (1, C('))- ^ e function z(-) is simply called a fluid model 
solution for (A, #, Co)- 

Note that 5(0, i) may be equal to +oo if Co = 0, i.e. z(0) = 0. Both right hand 
side terms in (|2.7|) are still well defined in this case, and the first term equals zero. 
The class of corner sets is defined as 

C = {[#, oo) x [y, oo) :i,j/G M+} U { [x, oo] x [y, oo] :i,i/£ R+} • 

The sets from the class C will be used to describe the evolution of fluid model 
solutions. Since each C = [x, oo] x [y, oo] £ C is characterized by the coordinates 

(x, y) of its corner, it is convenient to use the notation p{x, y) =' fx([x, oo] x [y, oo]) 
for any p, £ Mi. If zq > 0, for this class of subsets Equation (|2.7(l can then be 
rewritten as follows: for each x, y > and t > 0, 

(2.8) C(*)(*»V) = *oP(£° >x + 5(0,t),D°>y + t) 

, / P (B > x + S(s, t),D>y + t- s) ds, 
Jo 



+ A 
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Since z(t) — C(*)(0, 0), the fluid model solution z(-) satisfies the following equation: 
for each t > 0, 

(2.9) z(t) = z P (B° > S(0, t);D° > i) + A / P (B > S(s, t);D>t- s) ds. 

Jo 

It will be proved that the fluid model defined above is the limit in distribution 

of the rescaled processes {Z (•) : r € 1Z} introduced in Section ETT1 The measure 

valued fluid model solution £(•) corresponds to the measure valued state descriptor 

Z(-), and the fluid model solution z(-) is the limit of the queue length process Z(-). 

The main result concerning the convergence of {Z (•) : r e TZ} is the following 

theorem. 

Theorem 2.2. The sequence {Z (•) : r S 72.} is tight and each weak limit point 
is almost surely a measure valued fluid model solution £(•) for the data (A, $, Co)- 
// m addition l|2.10[l holds, then Z (•) converges in distribution, as r — > oo ; to i/ie 
unique measure valued fluid model solution C(')- 

This theorem is proved in Section [S] and 

2.4. Some Properties of the fluid model. Despite the quite abstract setting of 
this paper (measure valued processes), some concrete and explicit results concerning 
the fluid model of the queue can be obtained. Let (A, ■&, Co) satisfy the assumptions 
of Section IO 

The first result establishes the uniqueness of fluid model solutions under a Lip- 
schitz condition on the initial condition Co- 

Theorem 2.3. Suppose there exists a finite constant L such that 

(2-10) (o(Ax[y,y'})<L\y'-y\ 

for all Borel sets A o/l and all y' > y > . Then 1)2. 8JI and (|2.9(l have a unique 
solution. 

The second theorem analyzes the equilibrium of the fluid model, i.e. the behavior 
at infinity of the solution of Equation l|2.9l) . 

Theorem 2.4. Suppose that AE [-Bl.m =00 i1 < 1 and E [mh{B, D}] < oo. Then 
any solution t — » z(t) of l|2.9|l converges at infinity to the unique positive solution 
Zoo of the fixed point equation 

Zoo = AE [minjZooB, D}] . 

Moreover, any solution (C(i)) of Equation (|2.8|l converges to the measure Coo S Mi, 
defined by 

(2.11) C, 00 (x,y) = \[ P[B>x+—,D>y + t] dt 

Jq \ Zoo / 

= E [mm{zoo(B - x)+ , (D - y)+}] , 

for x,y > 0. 

These theorems are proved in Section|3J The simple fixed point equation stated in 
this theorem is used to analyze the qualitative behavior of the queue, see Section^] 
Note that the expression of the distribution of points Coo describing the asymptotic 
behavior of this queue has a simple expression in terms of the solution Zoo of the 
fixed point equation. 
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3. Some Properties of the Fluid Model 

In this section some basic properties of fluid model solutions are derived. In 
what follows, let z(-) be an arbitrary fluid model solution, i.e. such that 

fity z{t) = z P (B° > 5(0, t); D° > t) + A / P (B > S{s, t);D>t-s) ds. 

Jo 

If zq > 0, define 

S{t) = inf {s : 5(0, s) > i}, 

since z(t) < At + Zq, 5(0, t) — > oo as t — > oo, implying that S(t) is well defined for 
all t. In addition, S(t) < oo for all t if zq > 0, which follows from property 10) and 
continuity of the fluid model solution z(-). 

Define z(t) = z(S(t)), then (z(t)) satisfies the equation 



(3.1) z{t) = z P (B° >t;D° > S(t) 

+ A J z(u)P (B>t-u;D> S(t) - S(u)) du. 

We next introduce the concept of a shifted fluid model solution: For to > define 
z(to, t) = z{to + t), and define S(to, u, v) — J l/z(to, r)dr. 

Property 3.1. If z(-) is a solution of I2.fty . then z(to, •) = z(t,Q + •) satisfies 

z(t ,t) = C(t )(S(to,0, t),t) + \ V{B> S(t , s,t);D>t-s)ds,t>0. 

Jo 



Proof. Note first that by definition, and since S(to, 0, t) = S(to, to + t), 

C(«o)(S(*o, 0, t), t) = z P (B° > S(t , 0, t) + S(0, t ),D° >t a +t) 

/•to 
+ A / P (B > S(t - s, t ; D>s + t) + S(t , 0, t)) ds 
Jo 

= z P (B° > 5(0, to + t); D° > t + t) 

/•to 
+ X P(B>S{to-s,t +t);D>s + t)ds. 
Jo 

= z P (B° > 5(0, t + t); D° > t + t) 

/■to 
+ A / P(B > S(s,t + t);D > t +t-s)ds. 
Jo 

We apply this expression as follows: 

z(t ,t) = z(t + t) 

= z Q P (B° > 5(0, t + t); D° > t + t) 

+ A/ P(B > S(s,t + t)-,D > t Q + t-s)ds 
Jo 



fto+t 

C(to)(5(t , 0, t), t) + A / P (S > S(s, t + t); D > t + t - s) ds. 
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The change of variables y = s — to and the identity S(to + y, to + t) = S(to,y, t) 
give the result. □ 

The next proposition shows that continuity of fluid model solutions is a conse- 
quence of properties 101 and (JuJ . 

Lemma 3.2. Let the distribution of(B,D°) be free of atoms. Then any solution 
z{i),t > to 12. ,9|) satisfying inf t>a zlt) > for all a > is continuous. 

Proof. The function t -»■ 5(0, t) is continuous and so is t -> P (B° > 5(0, t);D° > t) 
since the distribution of (B°,D°) has no atom. The first term of the right hand 
side of Equation (|2.9|) is a continuous function of t. Concerning the second term 
of Equation i|2.9|) . by monotonicity the integrand t — ► P (B > S(s, t); D > t — s) is 
continuous almost everywhere on R + , hence its integral is a continuous function of 
t by Lebesgue's Theorem. The lemma is proved. □ 

3.1. A Maximal Solution. An important monotonicity property of fluid model 
solutions is proved in this section. 

Proposition 3.3. If AE [-B1{d =00 }] < 1 an d E [min{_B, D}] < oo, then any fluid 
model solution is bounded. 

Proof. Note first that, since E [min{i?, D}} < oo, also E [min{i?, aD}} < oo for 
every a G [0, oo). Define \\z\\t — s uPo<u<t z ( u )- Note that \\z\\t < zo + Xt < co. Fix 
t and let u £ [0, t]. Since S(u,s) > (u — s)/\\z\\ t , it holds that 

z(u) < z + X P{D>u~s; \\z\\ t B > u - s) ds < z + AE [min{L>, ||z|| t S}] 
Jo 

which is finite since E [min{£), B}] < oo. By taking the supremum over u £ [0,t] 
and by dividing both sides by ||z||t one obtains the relation 

l<zo/\\z\\ t + XE[mm{D/\\z\\ t ,B}}, 

If \\z\\t —>■ oo then, by monotone convergence, one gets the inequality 

1 < AE [Bl {D=oc} ] , 

which contradicts the assumption AE [-B1{d =0 o}] < 1- We conclude that \\z\\ t 
converges to some finite constant M which implies the assertion. □ 

Proposition 3.4 (Maximal Fluid Solution). For zq > 0, there exists a fluid model 
solution z*(-) starting form zo which is maximal, i.e. for any fluid model solution 
z(-) such that z(0) = zo, the relation z(t) < z*(t) holds for all t > 0. 

Proof. To define z*(-) we first define a sequence of functions z n (-), n > 0, by z°(i) = 

z + Xt, S n (u,v) = J2(l/z n (r))dr and 



l (t) = z P (B° > 5"(0, t);D a >t)+X f P (B > 5" 

Jo 



(t- s,t);D > s)ds. 



We show that z n+1 {t) < z n (t) by induction. The inequality z l (t) < z°(t) is trivial. 
Suppose now that z n (t) < z" _1 (i). Then S n (u,v) > 5" _1 (u,w), and, using the 
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fact that tail probabilities are non-increasing, 

z n+1 (t) = z P (B° > S n {0,t);D° >t)+X f P (D > s;B > S n (t - s,t))ds 

Jo 

< z P (B° > S 71 - 1 ^, t): D° > t) + A / P (B > 5 n_1 (t - s, t); D > s) ds, 

Jo 

which equals z n (t). 

Since z n (t) is decreasing in n and non- negative for all n there exists a function 
z*(t) such that z*{t) — linin^oo z n (t). By the definition of z n (t), we see that z*(t) 
satisfies (|2.9|) . 

Furthermore, we have z(t) < z*(t) for any given fluid model solution z(-). This is 
true because z(t) < z (t), and using an inductive argument as above, z(t) < z n (t) 
for every n. Since we know that at least one fluid model solution exists, it follows 
that inf t>a z*(t) > for every a > 0. By Lemma [3.21 it follows that z*(t) > is 
continuous. We conclude that z*(-) is indeed a fluid model solution. □ 

3.2. Convergence of fluid model solutions. In this subsection we show the 
convergence of fluid model solutions to a non-trivial constant Zqo as t — > oo. 

Proposition 3.5. If XE [B^{D=oo}] <1, E [min{S, D}] <oo and p>l, the equation 

(3.2) Zoc = AE [minors, D}] , 

has a unique solution in (0, oo). 

Proof. The function / : a — > AE [min{i3, aD}] is non-decreasing and concave on 
[0, +oo), note that f(a) = AE [min{B, aD}l{ D<+oa y] + AE [B1{_d=+oo}] for a > 0, 
therefore /(0+) = AE [Bl{ C=+00 }] <1 and f(a) converges to AE [B] >1 as a goes 
to infinity. By continuity of /, there exists do, 0<ao< + oo such that f(a) = 1. The 
concavity and the monotonicity imply that such a ao is unique, otherwise / should 
be constant equal to 1 after ao, but this is impossible since / converges to AE [B] >1 
at infinity. The quantity 1/ao is then the unique solution of Equation (|3.2I) . □ 

We are now ready to present the main result of this subsection, concerning the 
asymptotic behavior of any fluid model solution (z(t)) as t goes to infinity. 

Theorem 3.6. If z(-) is a fluid model solution, under the conditions 

AE [S1{d=oo}] < 1, E [min{B, D}} < oo and p > 1, 

then z(i)^Zoo as t—* + oo, where z^ the unique positive solution of the fixed point 
equation IS.ty) . 

Proof. It suffices to show z = limsup^^ z*(t) < Zoo and z = liminft^oo z*{t) > 
Zoo. We start with the former. We know that z < oo from Proposition ^. 31 For any 
e > there exists a t £ such that z(t) < z + e. We see that, for t > t £ , 

z(t) <z a P{B° >S(0,t);D°>t) 

rtts ft—t e 

+ X P(B>S(s,t);D>t-s)ds + X P (min{L>, (z + e)B] > s) ds. 

Jo Jo 

Taking the limsup on both sides, and noting that S(s,t) — > oo for any s > 0, we 

obtain that 

z < E [min{D, (z + s)B}\ . 
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The result is valid for every e > 0. By letting e [ 0, we obtain that z < Zoq. The 
lower bound follows by a similar argument after first noting that z > since z{-) is 
a fluid model solution. □ 

3.3. Uniqueness of fluid model solutions under non-zero conditions. The 

uniqueness of fluid model solutions is, in general, difficult to determine. If one looks 
at the time-changed version 1)3.1(1 . and take A = 0, one gets an ODE. Uniqueness 
of solutions to such ODE's can usually only be established by reducing it to some 
special case or to assume some kind of Lipschitz condition. If D = oo, then (|3.1() 
reduces to a renewal equation for which uniqueness is known to hold. Unfortunately, 
this reduction is not possible in general, which lead us to use a Lipschitz condition 
on the distribution function of {B , D ). It is not necessary to assume regularity 
conditions on the distribution of (B, D). We shall give a direct proof of uniqueness; 
for related results in the functional analysis literature, we refer to Chapter 2 of Hale 
& Verduyn Lunel |15|. 

Theorem 3.7. Suppose zq > and Fo(x,y) = P [B° > x: D Q > y) is Lipschitz 
continuous in Y , i.e. there is a constant L such that for any x,y,y', 

\F (x,y)-F (x,y')\<L\y-y'\. 

Then 12. 9\) has a unique solution. 

Defining ((t)(u,v) = <^(S(t))(u,v), it can easily be shown that 

(3.3) ((t)(u, v) = z a P (B >u + t; D° > v + S(tj) 

+ X J z(s)P (B>u+(t-s);D>v + S{t) - S{s)) ds. 

It is clear that for any u, v, ((t)(u, v), t > 0, is completely determined by z(t), t > 
and the initial measure. Thus, uniqueness of z(t) on an interval A carries over to 
uniqueness of £(i)(u, v) on A. 

The idea of the proof is simple: we take a suitable constant a > and prove first 
that uniqueness holds for z{t) for [0, a}. As discussed above, uniqueness carries over 
to CW( U > v ) f° r t £ [0, a]. Using this and the shifted fluid model equation given by 
Property 5.1, we prove uniqueness for z(t) on the interval [a, 2a], and so forth. This 
iterative procedure works if the measure £(£)(«, v) is Lipschitz for any < t < T. 
This is the content of the following lemma. 

Lemma 3.8. For any x, y, y' and for any t we have 

\((t)(x, y) - at)(x,y')\ < (z L + X)\y - y'\. 

Proof. We may take y, y' such that y < y' . From (|3.3|l we obtain 

\at)(x,y)~C(t)(x,y')\<z L\y'-y\ 

+ X J z(s)P (S(t) ~ S{s) +y<D< S{t) - S{s) + y'\ ds. 

Noting that z(s)ds — dS(s) we can rewrite this into 

rS(t) 

(3.4) z L\y'-y\+X P (r + y < D < r + y')ds 

Jo 
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Noting that, for any 5 > 

A / P (y < D < y + 5) dy < S, 
Jo 

we see that (|3.4|) can be upper bounded by (zqL + X)\y — y'\. □ 

Proof of Theorem \3.7\ By the one-to-one correspondence between solutions of (12.91) 
and (|3.1[1 . it suffices to show that 1|3.1|) has a unique solution. Define a — l/(2(zoL+ 
4A)). We first show that 1)3. Q has a unique solution on the interval [0, a]. For that, 
suppose that there exist two different solutions z(t), < t < a and h(t), < t < a. 
Set e = sup < t < a \z(t) -h(t)\. 
Note that for any < a < t < a, 

\H{t) - H{s) - (S(t) - S(s))\ < ea. 

Using (13 . 1|) for both z and h, together with the Lipschitz assumption, we obtain, 
after some simple estimates, 

\z{t) -h(t)\<z P (b° >t;D° > S(tj) - P (B° > t: D° > H(tj) 

+X J \z(s)P (b > t-s; D > S(t)-S(sj\ -h(s)P {B > t-s; D > H(t)-H(s)) \ds. 
The first term is bounded by z^Lae. The second term is bounded by 

A f \z(s)-h{s)\ds + 
Jo 

A / z(s) p(B>t-s;D> S(t) - S(sj) -P(B>t-s:D> H{t) - H{s)) ds. 

Call these terms Ila and lib. We have Ila < Xea. To bound lib, we use the 
bound 



P ( B > t - s; D > S(t) - 5(a) J - P (B > t - s; D > H(t) - H(s)) 

< P (S(t) - S(s) -ea<D< S(t) - S(s) + ea 



to obtain (after a change of variable r = S(s)) 

fS(t) 

IIb< X P (r - ea < D < r + ea) da < X2ea. 

Jo 

Putting everything together, we see that for t G [0,a], 

\z(t) - h(t)\ < z Lae + 3Xea < e/2, 

which implies that e = 0, i.e. that z(t) and h(t) coincide on [0, a]. Hence Equa- 
tion 1)3. ip has a unique solution in the interval [0, a]. 

Suppose now that (|3.1fl has a unique solution on [0, ka] for some k > 1, and 
consider the equation 

(3.5) z{ka, t) = Q{ka){t, S{ka, t)) 

r-l 



X / z(ka, a)P (B>t-s;D> S(ka, t) - S(ka, a)) ds. 
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We now show that this equation has a unique solution on [0, a], implying that there 
exist a unique solution of (|3.1(l on the interval [0, (k + l)a]. Suppose z(ka,t) and 
hit) both satisfy i|3.5fl . and set e = sup te r a i \z(ka, t) — h{t)\. As before, we have 

\(S(ka,t) - S(ka, s) - (H(t) - H(s))\ < (t-s)e <ae. 

Using this, we get as before (using now the Lemma for the first term) that 

\z(ka, t) - h(t)\ < {z L + X)ae + 3\ae < e/2, 

which implies that e = 0, and that uniqueness of solutions of l|3. ljl holds on the 
interval [0, (fc + l)a]. Iterating this argument completes uniqueness for all t. □ 

3.4. Uniqueness starting from zero. The result in this subsection can be seen 
as an extension of a result of |28| . who considered the case PS queue without 
impatience, i.e. with D = +oo. 

Theorem 3.9. Let e > 0. Suppose that (B,D) and a non-increasing function 
F e [x,y), with < F e {x,y) < Xe are such that 

z e {t) =F e (S s {Q,t),t) + X I P(B> S e (t-s,t);D>t-s)ds 

Jo 

has a unique solution z e {t) satisfying inf t>a z £ (t) > for a > 0. Then z £ (t) —* Z^(t) 
as e 10. 

Proof. As in the construction of the maximal fluid solution, z e {-) can be defined as 
the pointwise limit lim rwoo z™(-) with z™(-) recursively defined by z® = e + Ai and 

[ (t) =F e (S?{0,t),t) + \ f P(Z)>s;S>S , ; i (t-s,t))ds. 
Jo 

From this construction it can be easily shown that z™(t) is decreasing in n, and 
that z™(t) > Zq (t). Since also z*(t) < z 7 e l (t), we see that 

lim sup z e (£) < limsupz™(t) = Zq(*)' 

e|0 e|0 

Since this holds for any n, and Zg (t) — > z${t), we can let n — ► oo to obtain 

limsupz e (i) < Zg(i). 

To prove the other bound, we observe by induction and the properties of F £ that 
z o (t) < z s(t) for every n > 0. Consequently, 

z*(t) = lim z l (t) < lim sup z^(t) = z € {t). 

n >oo n — >OD 

We conclude that z e (t) > zo(i) for every e > 0, which implies the lower limit and 
the convergence z £ (t) — > z(t). D 

Uniqueness of fluid model solutions starting from is now a simple corollary. 

Corollary 3.10. Suppose that zq = 0. Then J^S.l)) has a unique solution. 

Proof. Let z(-) be a fluid model solution. Define z e (t) = z(t + e). Given z(s), < 
s < e, z e (-) satisfies the equation 

z e (t) =F e (S e (0,t),t) + A J P(D>s;B> S e (t - s,t))ds. 

Jo 



z" +1 ' 
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Here (with obvious notation) 

FJx,y)= / P\D>s + y;B>x+ / — -du ) ds. 

Jo \ Je-s z ( u ) J 

We see that F £ is globally Lipschitz in the second coordinate (with Lipschitz con- 
stant 1). Consequently, the above equation has a unique fluid model solution 
in terms of F e so that z e (-) is uniquely determined by z(t),0 < t < e. Since 
F e (x,y) < Ae, we see from the previous theorem that z e (t) — > z^t). But a l so 
z £ (t) = z(t + e) — > z(t), since z(t) is continuous. We conclude that z(t) = Zq(€), 
which implies uniqueness. □ 

— 2 
3.5. Analysis of the measure- valued fluid model. For any Borel set F of R , , 

the measure valued function £(•) satisfies the equation 

(3.6) t(t)(F)=to(F+(S(0,t),t)) + \ J d{F+{S{s,t),t~s))ds. 

Jo 

The properties, which are analogues of properties of z(-), are gathered in the fol- 
lowing theorem: 

Theorem 3.11. Let £(•) be a solution of fff.6)) . 

(i) Suppose p > 1,E [min{_B, Z?}] < oo and XE[Bl^ D=OQ y] < 1. ^4s i — > oo, 
£(£) converges to the limiting measure Coo defined by 



Co([x,oo] x [y, oo]) = A 



/■OO 

/ P (B - s >*/««>, £- tf >*) d«, 

Jo 



where z^ is the unique solution of the fixed point equation (|3.2|) . 
(ii) If Condition \2.1U\) of Theorem \2. !A holds. then Equation (|3.6|) has a unique 
solution. 

Proof. We know that z(t) — * z x as t — > oo. Consequently, S(t — to,t) —* to/ Zoo for 
every to > 0. Write for any Borel set F, 

C(t)(F) = C(0)(F + (5(0, £), t))+xf° *(F + (S(t -s,t), s))ds 

Jo 

+ X I tf(F+{S{t-s,t),s))ds. 
Jt 

Number the three terms on the right hand side as /, 77, III. By shifting time if 

necessary, we may assume that z(0) > 0. The first term converges to 0. Since 

z(0) > 0, there exists an r\ > such that z(t) > r\ for all t > 0. This implies that 

S(t — s,t) < s/rj. Consequently, since F C R 2 , 



III < A / P(B > s/m D > s)ds. 

Jtn 



'III 

From this bound, it follows that III — > as to - * oo. Since i? only has countably 
many discontinuities, and S'(t — s, £) — > s/z on [0, to] we have that 

II-»A / ^( J F + (s/z,s))ds. 
Jo 

Taking £ — > oo and then to ^ oo yield the first statement of the theorem. 
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To prove the second statement, note that Ij3.6(l has a unique solution for F = R 2 , 
which uniquely determines S(s, t) for all s, t with s < t. Since C is completely 
determined by Co: an d S(s,t), uniqueness follows. □ 

4. Applications 

In this section we analyze a number of quantitative properties of the fluid model 
equation l|2.9[l . In particular, we investigate the fixed point equation 

(4.1) zoo = AE [min{zooB, D}] . 

We treat a number of examples which allow for explicit computations, and also 
obtain a number of stochastic ordering results. In addition, we investigate the 
time-dependent behavior of z(t) for exponentially distributed lead times. 

We first give a heuristic interpretation of Equation (|4.1|1 : Let Z r denote the 
steady-state number of customers in the system. Furthermore, let V r (B) be the 
sojourn time of a customer if the customer never reneges. Then the actual sojourn 
time is given by mm{V r (B), D r }, and from Little's law we get 

(4.2) E [Z r ] = AE [min{V r (B), D r }] . 

Divide both sides of l|4.2[l by r. Since we observe the system in steady state at 
time 0, the number of customers hardly changes and by the snapshot principle 
we conclude that V r — Z r B + o(r). Furthermore, we have D r = Dr. Noting that 
Z r /r — » z then gives 14.1J1 after dividing both sides of Ij4.2[l by r and letting r — > oo. 
Apart from the mean queue length z, we are also interested in the long term 
fraction of customers that leave the system successfully. Denote this fraction by 
P s . It is clear that 

P s = P(D> Zoo B) . 

The following remarkable property, which simply follows from the fixed-point 
equation (|4.1|l , shows that the performance of the system does not depend on the 
average of D. 

Property 4.1. Consider two systems numbered by 1 and 2 such that (B2,D2) = 
(Bi,aDi) for some a > 0, and such that Ai = A 2 . Then (with obvious notation) 
we have 

Z2,oc = 0,Zl tOO , P s ,2 = Ps,l- 

We now proceed by analyzing a number of special cases. In Section 4.1, we 
assume a strong form of dependence. Section 4.2 assumes that B and D are inde- 
pendent. We give a remarkably simple expression for z(t) in the case that D has 
an exponential distribution. Finally, Section 4.3 considers an example which can 
be used as a flow level model for the integration of elastic and streaming traffic. 

4.1. Completely dependent lead times. Consider first the case where D — QB, 
with > (independent of B) reflecting the average service rate expected by a 
customer. In this case, the performance measures can be determined from the 
equations (recall that p = aE [B] > 1) 

Zoo = pE [min{6, Zoo}] , P s = P (@ > ^oo) • 

Some specific examples: 
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© single-valued. If we assume that = 8, then Zoo = pmm{8, Zoo}, which 
implies that Zoo = p8 since p > 1. From this, it follows that all customers 
leave the system impatiently: P s = P (6 > p8) = 0. Observe that when 
a customer leaves the system, a fraction 1/ p of his service time has been 
processed. 

two-valued. From the previous example, it is clear that the system can 
only get some work done if some customers are more patient than others. 
In this example we assume that equals 8\ with probability p and 82 with 
probability 1 — p. Take 82 > 9\. Equation (|4.1|) now simplifies to 

Zoo = ppmin{z oc ,8i} + p(l - p)min{z 00 ,8 2 }. 

From this equation and the properties 82 > 9i, p > 1 it follows that Zoo > 
6\. Furthermore, z x > 82 holds if and only if the equation 

Zoo = PP(>1 +p(l -p)Zoo 

has a non-negative solution, which is the case if and only if p(l — p) < 1 
(i.e. when the most patient customers cannot saturate the system alone). 
In this case we have 

"""' <8 2 . 



- l-p(l-p) 

If the last inequality is not valid or if p(l — p) > 1 we must have Zoo > 82 
which implies 

Zoo = PpOl +P(1 ~P)(>2- 

From the above we can conclude that P 8 = iff (1 — p(l —p))&2 < pp81. If 
the reverse inequality holds then all customers of type 2 are being served 
successfully, i.e. P s = (1 — p). 
- exponentially distributed. Assume w.l.o.g. that the mean of equals 1. 
In this case z^ can be determined from the equation z,^ = p{\ — c~ Zoc ) 
and P s = e~ 2 =° = 1 - z^jp. 

Since P s does not depend on the mean of 0, and since the worst-case property of 
the case of constant 0, it seems natural to conjecture that the system performance 
is positively related to the variability of 0. Thus it seems worthwhile to look for 

cvx 

ordering relations for P s if ©1 > ©2. If E [©i] = E[@2], this is equivalent to 
E [min{x, ©i}] < E [min{a;, © 2 }] for all x > 0. 

cvx 

Thus, if ©i > ©2, it follows that Z2 i0 o > 21,00 i.e. less variability in reneg- 
ing behavior implies a lower service rate. To prove that also P (©i > zi i00 ) > 
P (©2 > Z2,oo) seems hard without imposing further assumptions. 

4.2. Independent lead times. In this case we can write (|4.1|) as 

/•OO 

A/ P (B > u) P {D > Zoou) dw = 1. 

which, in case E [B] < 00, is equivalent to P (D > z^B*) = 1/p, with B* a random 
variable with density P (B > x) /E [B] . 

Recall that P s = P (D > z^B). Consequently, if B is exponentially distributed, 
we have the insensitivity (w.r.t. the distribution of D) result P s — 1/p. The in- 
equality P s < 1/p holds if B* is stochastically dominated by B, and P s > 1/p vice 
versa. Since B* being stochastically dominated by B is related to a low variability 
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of B, we see again that more variability (this time in the service times) leads to a 
better system performance (i.e. higher P s ). 

Exponential reneging 

If we assume that D has an exponential distribution (and B a general distribution) , 

we see that Zrx, is the solution of 

(4.3) p[3*{z 00 v))=l, 

with P*(s) = E[e- sB *]. 

In addition, we have the following remarkable expression for the complete fluid 
limit z(t),t > 0, if z = 0: 

Proposition 4.2. Suppose P (D > t) = e~ ut , that B is independent of D and that 
zq = 0. Then the unique solution of \2.ty is given by 

(4.4) z(t)=z 00 (l-e- ut ), 

with z the solution of Equation \4-?fy - 

Proof. Recall that Equation l|2.9|) has a unique solution. We show that i|4.4|) is 
indeed the solution of 1)2. 9J1 by verification. We thus compute the right hand side 
of (E2J) writing 2; (it) = Zoo (l - e~ vu ). 
Observe that 



, f -}-du = log(e^ - 1) - log(c- - 1). 



Consequently 

ft 



A / P (D > t - s)P ( B > / (l/z{u))du\ ds 

= -e" 1 " / P ( Zoo B > log(e^ - 1) - log(e" s - 1)) de vs 
v Jo 

= -e' vt / e—P { Zoo vB > log(e ,/i - 1) + Zoo ) dv 

v J-logie" 1 -!) 



■ log(e" 

-e- vt (e vt - 1) / P (ZoovB > v) e- v dv 
v Ja 



= Zoc (l - e- vt ) P P*( ZooV ) = Zoo (l - e"^). 
Which shows that Zoo(l — e~ vt ) satisfies (12. 9|) . 



D 



4.3. TCP-friendly traffic. Assume that there exist independent random vari- 
ables B\ and D\ with finite means such that 

(B,D) = (I?x, 00) with probability p, 

= (co, Di) with probability 1 — p. 

When we view PS as a way of modeling TCP, this example models the integration 
of elastic (TCP) traffic and TCP friendly UDP traffic; see Key et al. QH for a related 
model. The latter type of traffic is using the system for a certain amount of time, 
regardless of the level of congestion. 

The fixed point equation (|4.1|l for q specializes to 

zoo = ApE [zooBi] + A(l - p)E [£>i] , 
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Consequently, if the stability condition XpE [Bi] is satisfied, we see that 

A(l-p)E[£>i] 
Z °° l-ApE[Bi] ' 

5. Tightness 

In this section we prove the first part of Theorem 12. 21 that is, we show that the 
sequence of processes {Z (■), r G 1Z} is tight in D([0, oo), Mi). The main results in 
this section implying this property are the compact containment Lemma |5. 21 and 
an oscillation inequality in Lemma 15.61 To prove these results, a number of further 
lemmas are developed. Section 15.11 derives a Glivenko-Cantelli theorem for the 
stochastic primitives. Section 15 . 21 introduces a fluid scaled version of the dynamic 
equation for Z (•). The compact containment property is derived in Section 15.31 
Section 15.41 serves as a preparation for the oscillation bound. In particular, it 
is shown that Z (t) charges arbitrarily small mass to thin L-shaped sets. The 
oscillation bound is then shown in Section 15.51 

Throughout this section, it is assumed that the assumptions of Section l2~2"l hold. 

5.1. A Glivenko-Cantelli theorem. An important preliminary result is the fol- 
lowing functional Glivenko-Cantelli theorem for the stochastic primitives. It will 
be convenient to consider them together as a single, measure valued arrival process. 
For r G 1Z and t > s > 0, define the fluid scaled measure valued arrival process by 



£(*) = - 2^ <5(Br,13J>-i), 



4=1 

and define the fluid scaled increment 

(5.i) T{s,t)=T{t)-T{s). 

Note that C (•) is a random element of D([0, oo),Mi) and, for each t > s > 0, 
C (s,t) is a random element of Mi. 

To state and prove the result, we first introduce some notions from empirical pro- 

2 

cess theory. Our primary reference is 32 . A collection C of subsets of M + shatters 

2 

an n-point subset {xi, . . . , x n } C R + if the collection {C (~l {xi, . . . , x n } '■ C G C} 
has cardinality 2™. In this case, say that C picks out all subsets of {x\ 1 . . . ,x n }. 
The Vapnik-Cervonenkis index (VC-index) of C is 

Vc = min{n : C shatters no n-point subset}, 

where the minimum of the empty set equals infinity. The collection C is a Vapnik- 
Cervonenkis class (VC-class) if it has finite VC-index. 

VC-classes satisfy a useful entropy bound. Let Q denote the set of Borel proba- 

2 

bility measures on R + and, for Q G Q, let ||/||q = (|/|,Q) denote the Li(Q)-norm 

2 

of a Borel measurable function / : M + — » R. For e > 0, the L\{Q) e-ball around 
/ is the set of Borel functions {g : \\f — g\\Q < e}. For a family of functions Y, 
the (e,ii(Q))-covering number JV(e, ~f,Li(Q)) is the smallest number of L\(Q) 
e-balls needed to cover "V . If C is a VC-class, then for all e > 0, the family 
f = {\ C ■■ C G C} satisfies 

(5.2) suplog7V(e,r ) L 1 (Q))<oo; 

QeQ 
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see Theorem 2.6.4 in [52"] . 

Recall the collection of corner sets C defined in Section l2~3l 

C = {[x, oo) x [y, oo) : x, y G M.+ } U { [x, oo] x [y, oo] : x,y £ K+} • 

2 

Note that for any 3-point subset {a;i,X2,a;3} C R + , it is impossible for C to pick 
out all three 2-point subsets of {x\,X2 1 x^\. Since C shatters no 3-point subset, it 
has VC- index bounded above by 3. Thus, C is a VC-class and y = {lc ■ C E C} 
satisfies (j5.2|) . 

Define an envelope function for Y as follows. Let 7r : R + —> R+ be the map 
n(x,y) = max{x, y}. Since ir is continuous, l|2.5|) and the Skorohod representation 
theorem imply the existence of K + -valued random variables X r ~ ?? r o ir^ 1 and 
X <~ i? o 7T _1 such that X r — > X almost surely. Thus, there exists an R + -valued 
random variable Y such that 



(5.3) 



Y = supX r , almost surely. 

r£K 



Let \x be the law of V on R + . Since L,2(p) contains continuous unbounded functions, 
there exists a continuous, unbounded function ip : R + — > M + that is increasing on 
[0, oo), satisfies ip > 1, and such that (ip 2 , fi) < oo. This implies that 



(5.4) 



((tp o tt) 2 ,^) = E[i/;p0 2 ] < E[i/;(r) 2 ] < oo. 



Let F = %p o 7r, and note that lc < -F 1 for all C G C. That is, F is an envelope 
function for ^. Finally, define f = f U {F}. 

Lemma 5.1. Let T > 0. Then as r — > oo. 



(5.5) 



sup sup 
/eT0<s<t<T 



(/,r( s ,i))-A''(t- s )(/,o 



Proof. Let e > 0. By l|5.1|l . it suffices to show that 



limsupP r sup sup 
r->oo V/eT*e[o,T] 



(f,C r (t))-\ r t(f,0 r ) 



> e\ <e. 



Note that the above event is measurable for each r because it can be rewritten 
using the suprema over rational t, and / = lc with C having rational or infinite 
corner coordinates x and y. Since (/, C (£)) and \ r t(f,d r ) are nondecreasing in t 
for each fixed / G "f, it suffices to show that for each fixed t G [0, T], 



limsupP r sup 



(/,£ (t))-A r i</,0 p ) 



> e < e. 



Since 



(/, C r (t)) - X r t(f, #") - (/, #■) (V(t) - A r t) + -ff{t) 



T{t) 



(f,$ r ) 
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(with the convention that division by zero equals zero), it suffices to show the two 
bounds 



limsupP r sup 

r-,00 \fer 



(5.6) limsupP r sup 



(fJ r )(E r (t)-X r 

T* n f(fXV)) 



E{t) 



E{t) 



-</,^> 



e \ e 
>2 % 



e \ e 



The first equation follows from assumption 1)2. 4|l and by observing that 
(5.7) sup sup (/, ti r ) < sup(F, tf r ) = sup E[<0(X r )] < E[V>(F)] < 00, 

reTZ f^f rETl rETl 

which follows from l|5.3(l and l|5.4|) . To show l|5.6[) , it suffices to verify three assump- 
tions of Theorem 2.8.1 in [22]- Observe that for each n £ N and (ei, . . . , e n ) £ R™, 
the function 

n 

(xi,...,x n ) -> sup y^ ej/fo) 
/er i= i 

is measurable on the completion of (M. + ,3S, d r ) n , for each re R. Thus, Y is a 
d r -measurable class for each r E 1Z; sec Definition 2.3.3 in [22]- Moreover, ^ is 
uniformly bounded above by the envelope function F, and 

lim sup(Fl {F>M} ,& r ) =0, 

Af-KX) rG 7j 

by Markov's inequality, (|5.3|l . and (|5.4|l . Lastly, ^ satisfies the finite entropy bound 
JS3 because iV(e,r,Li(<9)) < i\T(e,r,Li(Q)) + 1 and C is a VC-class. The 
previous three observations imply that the assumptions of Theorem 2.8.1 in |32| 
are satisfied. Consequently, "V is Glivenko-Cantelli, uniformly in r. That is, for 
every S > 0, there exists an n$ such that n > n$ implies 



(5.8) 



supP' 

rETl 



I sup S up -]T/(5[, DJ-r" 1 ) - (fJ r ) >S)<5. 



Choose S = min{e/2, e/(4AT)}. The left side of i|5.6[) is bounded above by 



lim sup P 



(V(i) 



> 2AT) +limsupP r sup 

r^oo \feT 



(fx r (t)) 



E (t) 



(fJ r ) 



> 



4AT 



The first term equals zero by (|2.4I) . For the second term, rewrite 



(fX®) 



1 



E*-(rt) 



E r (rt) 



EH) - 

and bound each probability in the second term by 



Y, /(BlDlr- 1 ), 



(5.9) P r (E r (rt) < n s ) 



( sup sup - J2 m, DVr- 1 ) - (/, #■) > -i- j 



By (I2.4|l . the first term in l|5.9|l converges to zero as r — > 00. By 1)5. 8|l . the second 
term is bounded above by 5 < e/2, uniformly in r e K. This implies l|5.6[) . D 
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5.2. Fluid scaled dynamic equation. Using (|2.3(l . it is easy to see that the 
fluid scaled state descriptor of the rth model satisfies the following equation almost 
surely: for each Borel set A £ 38, and all t, h > 0, 

(5.10) Z r (t + h){A) = Z^(t) U + (S r (t, t + h), h)) 

l£(BT(t + fc),:z£(t + /»)). 



rE (t+h) 

E 

i=rE r {t) + l 



Subsequent proofs use estimates obtained from this equation. Two estimates 
result from bounding the summands in (|5.1U|) by 1 and optionally bounding the 
first term on the right side by its total mass; for each A £ 38 and t, h > 0, 

Z r (t + h)(A) <Z r (t)(A+(S r {t,t + h),h)) +t(t,t + h)(R 2 + ) 

( 5 - U ) -, )-2X -r -2 

<z (*)(»;) +c (t,t + h)(R + ). 

Two more estimates follow from (|5.1U|I by simply ignoring any arrivals; for each 
A £38 and t, h > 0, 

(5.12) Z r (t) (A + (S r (t, t + h), hj\ < Z r (t + h)(A) < Z r (t + h) (r 2 + \ . 

5.3. Compact containment. This section establishes the compact containment 
property needed to prove tightness. 

Lemma 5.2. Let T > and r\ > 0. There exists a compact set K C Mi such that 

(5.13) liminf P r (if \t) £ K for all t £ [0,T]) > 1 - r). 



2 

Proof. A set K c Mi is relatively compact if sup£ GK £(M + ) < oo, and if there 

2 2 

exists a sequence of nested compact sets K n C M + such that UneN ^ n = ^+ an d 

lim sup £(#£)=(), 

where K^ denotes the complement of K n ; see [TJ\, Theorem A 7.5. Consider the 

2 

nested sequence of compact sets in R + given by 

K n = ([0,n] x [0,n])U([0,n] x {oo}) U ({oo} x [0,n])U({oo} x {oo}), neff. 

By (|2.6|l . Z (0) — ► (o in distribution, and so the sequence {Z (0)} is tight. 
Thus, there is a compact set Ko C Mi, such that 

(5.14) lim infP r (if (0) S K ) > 1 - -. 

r — >oo 2 
2 

Let M = sup ?GK() £(R + ), an( i 1^ «n = su PfeK £(-^n) f° r eacn "E^. Since K 

2 

is compact, M < oo and there exists a sequence of nested compact sets J„ C R, 

2 

such that U„ eN J n = R + and lim^oo sup ?GKo £(J£) = 0. Since J n C i^fe(„) for 
each n £ N and sufficiently large fe(n) £ N, it follows that a n — ► as n — * oo. 

Recall the definition from Section I5~T1 of the envelope function F — tp on for the 

family T. By (23} and JE3, the constant M = sup,, eTC (\ r T(F, d r ) + l\ is finite. 
Let K be the closure of the set 

U e Mi : £(!+) < M + M and £(K%) < a n + V>(w) _1 M for all n £ n| . 
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Since a n + '0(n) _1 M — > as n — > oo, the set K is compact in Mi. 

For each r G 72, denote the event in (|5.14|) by fig an d define the event 



SV, 



J = {(F.ZTCT)) < X r T(Fj r ) + l} 



By (|5.14|) and Lemma Ol lim inf ,._>«, P r (r2o n fi£) > 1 — 77. Fix w G fig n QJ and 
£ G [0, T], and assume for the remainder of the proof that all random objects are 
evaluated at this u>. Then it suffices to show that Z (£) G K. 



By (EH, 



Z r (t)(R 2 + ) < Z r (0)(K+) + C r (t)(R 2 + ). 



Since £ r (i)(M + ) = (1,£ (*)) < (1,£ (T)) < (F,£ r (T)}, the definitions of fig, Of, 
and M imply that 

(5.15) Z r (t)(l+) <M +M. 
Fix n G N. By lf5~TU|l . 

z r (t)(^) =r(o) (iq + (s r (o, ;),*)) + E ^ (bT(*),sT(«)) ■ 

4=1 

The shape of the set TiT^ implies that 

#° + (5(0,t),t) C K° and 1+, (^(t),^(t)) < ^.(flf.DJV- 1 ), 
for i= l,...,rlf(t). Thus, 

Z r (i)(^ c J < Z r (0) (1Q) + <l*o,r(t)>. 
By definition of V', i*\ and by Markov's inequality, 1^-c < ip(n)~ 1 F. So 

Z>)(i^) < Z r (0) (7^) + ^(n)- 1 ^^^)}. 
Since (F,T '(*)) < (F,ZT(T)), the definitions of Og, fi^, and M imply that 

(5.16) ^"WTO < «n + ^W^M. 

Equations jETTS} and l|niini imply that z'ffjeK. □ 

5.4. Asymptotic regularity. The second and main step necessary to prove tight- 
ness is to bound the probability that the process Z (•) oscillates. Oscillations may 
result from sudden arrivals or departures of a large amount of mass. Sudden ar- 
rivals are controlled by the regularity of the arrival process. To show that sudden 
departures are unlikely as well, we show that Z (•) assigns arbitrarily small mass 
to the boundaries of the sets C G C. This is phrased in terms of re-enlargements of 
the boundaries of these sets (forming a collection of L-shaped sets) . For C G C and 
K > 0, let dc denote the boundary of C and let 

d c = i w e WL : inf ||te - z\\ < K 

2 

be the K-enlargement in K + of its boundary, where the infimum over the empty set 

2 

equals 00. (Note that dc, and therefore also d c , is empty for the corner sets K, 

and {00} x {00}. Note also that d c = ((x — k) + ,x + re) x {00} for a corner set of 

the form [x, 00] x {00} with x G [0, 00).) The following lemma establishes the result 

for the initial condition Z (0). 
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Lemma 5.3. For all e,rj > there exists a k > such that 

(5.17) liminf P r ( sup 3f(0)(dg) < e ) > 1 - 7). 

r ^°° \cec J 

Proof. Fix e,T) > and let Z[(0)(-) =I r (0)(-xl + ) andZ^(0)(-) = if(0)(R + x ■)• 
For each C & C and k > 0, 

<9£ C {[x, x + 2k] x 1+) U (1+ x [y, y + 2k]), 

for some (x, y) e R^ = [0, oo) x [0, oo). Thus, it suffices to show that, for i = 1, 2, 
there exists a k > such that 

(5.18) liminf P r [ sup %(0){\x,x + 2d) < e ) > 1 - -. 

r ^°° \xe[o,oo) J 

We prove the statement for j = 1; the proof is identical for i = 2. 

The projection (a;,y) h- > x is continuous, so (|2.6|l implies that 2^(0) converges 
in distribution to £o(" x R + ) as r — > oo. Since ( (- x R+) is free of atoms in [0, oo), 
there exists a k > such that 

(5.19) sup Co([x,x + 4k} xl + ) <e. 

a:6[0,oo) 

(If H5.19|l fails, it is easy to construct an atom of (q(- x R+).) Moreover, there exists 
a constant M such that 

(5.20) (o([M,oo) xR+) <e. 

Let N = \M/k] + 1, where \x] denotes the smallest integer n > x. For n = 
1, . . . , N — 1, define the set 7„ = [nn, (n + 4)k] and define In = [M, oo). Note that, 
for every x £ [0, oo) there is an n < N such that [x, x + 2k] C I n - To prove Ij5.18|l . 
it therefore suffices to show that 

(5.21) liminf P r ( maxlj(0)(/ n ) < e ) > 1 - -. 

r^oc yn<N J 2 

Let Mi(R_|_) denote the space of finite nonnegative Borel measures on R + , en- 
dowed with the weak topology. Let A = {£ 6 Mi(R+) : m&x n <N £,(I n ) < e}, and 
suppose that a sequence {£&} C Mi(R + ) satisfies £& — ► £ for some £ G A. Since 
the sets I n are closed, the Portmanteau theorem (adapted to finite measures) im- 
plies that 

lim sup £&(/„) < £(/«) < £, for all n < N. 

k — >oo 

Hence, £& € A for sufficiently large fc, which implies that A is open in Mi(R + ). 
Thus, a second application of the Portmanteau theorem yields 

liminf P r (z[(0) £ A) > P(Co(- x R+) £ A) = 1, 

r — >oo 

which implies l|5.21|l . D 

The regularity result is now shown for the entire state descriptor Z (■). 
Lemma 5.4. Let T > and e, r/ > 0. There exists a k > such that 

(5.22) liminf P r ( sup sup T{t){d£) < e ) > 1 - r). 

r ^°° \ cec te[o,T] I 
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Proof. By Lemmas l5.il 15.21 and 15. 31 there exists a compact K C Mi and a kq > 0, 
such that for all 6 > 0, the events 



Q{ = i supZ r (0)(^°) < £ 
.cec 



C- 2 



r?2 = \ su p su p 

lceco<s<t<T 



Z r (s,i)(C)-A r (i-s)i? r (C) 



<<5 



I = JZ 7 ^) 6 K for all t e [0, T)\ 



n5 = n;nfi5nn5, 

satisfy 

(5.23) liminfP r (£^) > 1 - 77. 

r — >oo 

Recall the compact sets i"T„ defined in the proof of Lemma l5~21 Since K is compact, 
there exists a finite M > 1 and an integer R < 00 such that 

(5.24) sup £(!+) < M, 

(5.25) SU p^(^)<£. 

CeK 2 

Let A* = sup reTC A r , which is finite by l|2.4[) . Fix 

h = e(8X*)-\ K = mm{n ,h(2M)- 1 } and 5 = e min{(8\ RM h,- 1 })- 1 , 2" 1 }. 

For r £ 1Z, let f2* denote the event in l|5.22|) . By (|5.23() . it suffices to show that 
f2p C fi*. Let w G 0[j be arbitrary; for the remainder of the proof, all random 
objects are evaluated at this u. 

Consider any r E 11, t e [0,T] and C e C. We must show that T (t)(d£) < e. 
Define the random time 

n =sup{s <t: (l,Z r (s)) =0}, 

if the supremum exists, and define t\ = otherwise. Let r = max{ri,i — RM}. 
We first show that 



(5.26) Z r (r)^ + (S r (r,t),t-r))< £ -. 

If t — 0, this follows from the definition of f2J because k < kq, because 

and because C is closed under positive translation. Suppose r = t± > 0. Then there 
is a sequence {r n }, with r„ f r, such that (1,-Z (T n )) = for all n. In this case, 
(I5.11|l and the definition of fl^ imply that, for all n, 

Z r (r) (dg. + (S r (T, t),t-r))< Z r (r n ) (l+) + T(t„,t) (l+) < X r (t - r n ) + S. 

Letting r n \ r yields 

Z r (r)^ + (S r (r,t),t-r))<5<^. 
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Suppose that r = t — RM. Since (1,2 (s)) > for all s £ (r, t], the definition of 
fig and (|5.24|) imply that 

ft 



It-RM 

Thus, by the definition of fig and 15. 2511 . 



S r (r,t) = f (l,2 r (s)y 1 ds>R. 

Jt-RM 



2\t) (ag + (S r (r,t),t - r)) < Z» (ff£) < |, 



which proves 1)5.26(1 . 

By (prrnf . 



(5.27) 2 r (t) (flg) = Z r (r) (c£ + (S r (r, t), * - r)) 

+ - E lfe(Bl(*).^(*))- 



=rE (t) + 1 



Let I denote the second right hand term in 1)5.27(1 . By 1(5.26(1 . it remains to show 
that I < e/2. Let N = \(t- rfh' 1 ] and, for each n = 0, . . . , N - 1, let £„ = r + nh 
and i" = rnm{t„+i, £}. Then, using the inequality 1q k (•, •) < lg* (■,■), 



JV-l r£ (t") 

n=Q i=rW'(t n )+l 



JV-l rJ3 (t") 

(5-28) '<E- E lag (^W^*) 



Consider n £ {0, . . . , N — 1} and i such that t/J> x G (£„, £"]. Observe that 

(5.29) ~S r {t n ,t) KS" \Ulr~ 1 ,t) <~S r {t n ,t). 
By definition, 

(5.30) loziKit), Tf&)) = la^^urr-^t-u^r-^BlDrr- 1 ). 
So, letting 

c~ = c+ (s r (t n , t) - k, t - r - k\ n i 

C„ = C + f S 1 (£„, t) + k, t — t„ + k J i i jr\ 2 . 
w = (7 n \ C n , 
it follows from 1(5.29(1 and 1(5.30(1 that 

(5.31) i a g(^(i),zr( f )) < i^B^r" 1 ). 

Conclude from (|F^5|) and (|53Tjl that 

N-l , rIT(t") 



■"C, 



nr H 



^(r(t n ,r)(c-)-r(t„,r)(c^ 



r _ 

n=0 i= r S r (t„)+l 



iV-1 



n=0 
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For all n < N, C~, C+ G C and t n — t n < h. So the definition of £l 2 implies that 

JV-l 

I < J2 (^h$ r (Cn) + 25) . 
By definition of N, and since t — t < RM, 

JV-l 

I <X*hJ2 ^ r (C„) + \RMh~ 11 [ 25. 
This implies, by choice of 5, that 

N-l 

(5.32) / < X*h Y, 3 r (Cn) + |. 

n=0 

If n€{0,...,iV-3}, then 

~s r (t n+1 ,t n+2 ) >hhr l >2k, 

because < (1,2 (s)} < M for all s G (r, i] and because h > k2M by definition. 
Thus, for all n e {0, . . . , JV - 3}, 

S T (t 7 \t) - K = S T (t n+ l,t n+2 ) + S r (t n+2l t) - K > S T (t n+2 ,t) + K. 

Hence, C~ C C^ +2 for all n G {0, . . . , N — 3}, and consequently, C n H C n +2 = 0- 
Thus, since $ r is a probability measure, 

L(AT-1)/2J L(W-2)/2j 

53 tf''(C 2n )and Y, ^(Csn+i) 

n=0 ra=0 

are both bounded above by one. Conclude from l|5.32|l that 

I <2X*h+^ 

which implies, by choice of h, that / < e/2. D 

5.5. Oscillation bound. This section establishes the second main ingredient for 
proving tightness of the state descriptors. As a metric on Mi, we use the Prohorov 
metric (adapted to finite measures). For fj,,v G Mi, define 

d\fi,u] =Mi e>0 : ^{A)<v(A e ) +e and i>(A)<(i(A 6 ) + e for all closed A E & 

Recall that A e = {w G R + : mi ze A \\z — w\\ < e} and that £3 denotes the Borel 

2 

subsets of K + . 

Definition 5.5. For each £(•) G D([0,oo),Mi) and each T > 6 > 0, define the 
modulus of continuity on [0,T] by 

W T («•).*)= SUp SUp d[C(* + ft),C(*)]. 

te[0,T-S] hG[0,6] 

Lemma 5.6. For all T > and e, r\ G (0, 1), i/iere exists S G (0, T) smc/i that 

(5.33) hminf P r (w T (z r (-), s) < e) > 1 - v. 
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Proof. As before, let A* = sup r(EK A r . For each k > 0, define 

L K = ([0, k] x R + ) U (1+ x [0, k]). 
By Lemmas 15. fl and 15.41 there exists n G (0, f) such that for all 8 6 (0,T), the 



events 



n\ = { sup Z r (t)(L K ) <-}, 
\te[o,T] 4 1 



n r 2 = 



sup C (t,t + 6)(R + ) <2\*S\ , 



n r = n[nn r 2 , 



satisfy 
(5.34) 



limmfP r (f^) >l-r). 



Fix S = K£ 2 (8max{A*, l})" 1 and let £11 be the event in l(5~33j) . By (|534) . it suffices 
to show that f2g C il* for each r. Fix r E 1Z and w e ilg; for the remainder of the 
proof all random objects are evaluated at this uo. Fix t e [0, T — 5], /i£ [0, S] and 
let A £ 3$ be closed. It suffices to show the two inequalities, 



(5.35) 
(5.36) 
To show 

(5.37) 



Z r (t)(A)<Z r (t + h)(A e )+e, 
Z r {t + h)(A) <Z r {t){A £ )+e. 
5j) , use the definition of Vt\ to write 

Z r (t)(A)<Z r (t)(L K ) + Z r (t)(AnLl) 
< £ -+T\t)(A^Ll). 



Let I = {s e [t,t+h] : (1,Z (s)} < e/4}. Suppose 7 = 0. Then (1,Z (s)) > e/4 
for all s G [£,£ + /&], which implies that 

rt+S —r 4(5 

(l,Z r (s))~ 1 ds + (5 < \-S < min{e, k}. 



(5.38) 



(5 (t,t + h),h) 



< 



it £ 

Consequently, (re, i/)e4nL^ implies (x, y) — (S (t, t + h),h) g A e , and so 
(5.39) AnL c K cA £ + (S r {t,t + h),h). 

Deduce from (|5.37() that 



Z (t)(A) <~ A +Z {t)[A' + (S {t,t + h),h) 



Apply (|5~T21) to get 
(5.40) 



Z r (t)(A)< £ -+Z r (t + h)(A"). 



Suppose I ^ and let r = inf I. Then (1,Z (r)) < e/4 by right continuity. Since 
(l,if (s)} > e/4 for all se [£,t), 



(5.41) 



(S (t,T),T-t) 



! (l,Z r {s))- 1 ds + 6<— + 5<K. 
t ' £ 



By (t5~37jl and l|OT)l . 



2 r (t)(A) < | +Z r (i)(L^) < | +Z r (t) (< + (S r (t,r),r-i) 
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Apply (f5~T5|> to get 

(5.42) Z r (t)(A)<|+r(r)(S+)<|. 

So H5.35[> follows because either (|5.4UI) or (|5.42[) holds. 



To show ()5.36[) , use (|5.11|l and the definitions of fi^ and 5 to obtain 

Z r {t + h)(A) <Z r (t)(A+(S r {t,t + h),h)') +T(t,t + h){W 2 + ) 

<Z r (t)(A + (S r (t,t + h),h)) + ^. 



(5.43) 



If J = 0, then (|5~35|) implies that A + (S" \t, t + h),h) C A e . So fB~4^)l yields 

Z r {t + K){A)<Z r {t){AZ)+ £ -. 
If I 7^ 0, then by (|5.11|l . the definition of f^ an d the choice of (5, 

Z r (t + h){A) <Z r (r)(I+)+Z r (r,t + / l )(I+)< |+2A*<5<|. 
In both cases, l|5.36[l holds. Conclude from 15.35fl and (|5.36l) that 



d 



<£. 



Z (t),Z (t + h) 
Since t € [0,T — 6] and ft, € [0,(5] were arbitrary, 

w T (z r (-),<5) <e, 
which implies that w G f2*. □ 

6. Limiting Fluid Equations 

This section contains the proof of Theorem 12.21 Tightness of the sequence 
{Z (•)} follows immediately from Lemmas 15.21 and 15.61 Since {Z (•)} is tight, 
there exists a subsequence {q} c 1Z and a process Z(-) in D ([0, 00), M) such that 
Z (■)=$■ Z(-) as q — ► 00. We must show that Z(-) is almost surely a measure valued 
fluid model solution for the data (A, $, Co)- This is accomplished by Lemmas 16.11 
and 16.21 and Theorem 16.31 below. Finally, if 1)2.10(1 holds, then a measure valued 
fluid model solution for (A, $, £o) is unique by Thcorem l2.3l In this case, the law of 
the limit point Z(-) is unique and so Z (•) =>■ Z(-) as r — > 00. 

Let Z(-) = (1,Z(-)} be the total mass process for Z(-), and let S(u,v) — 
lu ~zTsi ^ s ^ or a ll u > w > 0. To show that Z(-) is almost surely a measure valued 
fluid model solution, note first that Z(-) is almost surely continuous by Lemma 
15.61 Note also that, by l|2. 6J1 . Z(0) — (b almost surely. It remains to show that 
properties (Q) and (JuJ) of Definition ^ . II are satisfied almost surely by Z(-). The next 
result establishes (Q). 

Lemma 6.1. Almost surely, for all a > 0, 

(6.1) miZ(t)>0. 

t>a 

Proof. Take t > 0. Pick a constant a < t small enough such that the marginal 
distribution of D is continuous at a, take m < 00 such that the marginal distribution 
of _B is continuous at m, and such that AE [Bli D>a .B<m}\ > 1- By dominated 
convergence, 

lim A r E [i?[l{_DJ r -i> a; _BJ'< TO }] = AE [-Bl{£»a;B<m}] • 
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Compare the original system with an ordinary PS queues having arrival rate \ r am — 
A r P (D r > ra;B r < to) and service times BT am , which are distributed as B[ \ 
D\ > ra;B\ < to. Suppose that this PS queue is empty at time r(t — a), and let 
Z r (t) be the queue length in this PS queue at time ra. 

Observe that the number of arrivals in the modified PS queue between time 
r(t — a) and time rt is less than or equal to the number of arrivals in that interval 
in the original PS queue with impatience. Furthermore, if one of the jobs that 
arrived in the original PS queue after time r(t — a) departs before time rt, then this 
must also be the case in the modified PS queue, since that PS queue had a service 
rate which was at least as large as in the original PS queue. These considerations 
imply that Z r (rt) > Z r (rt). Since the modified queue is still overloaded, and no 
customer departed because of impatience, and the modified arrival process is still 
a renewal process, the evolution of the modified system between time r(t — a) and 
rt has the same law as that of an overloaded GI/GI/1 PS queue starting at 0, in 
the time interval [0, ra]. 

Since the service times in our modified system are bounded, the means converge. 
The assumptions in [2H| are therefore valid, and it follows that there exists a con- 
stant k a > such that linv^^ Z r (rt)/r = k a almost surely. Consequently, we have 
limnuv^oo Z (t) > m a almost surely, which implies the assertion. □ 

Before establishing property JnJ of Definition 12. II the following result is needed. 

Lemma 6.2. Almost surely, for all C G C and t > 0, 

(6.2) Z{t){d c ) = 0. 

Proof. Let T > 0. It suffices to show the statement for all t € [0,T]. Let {rj n } C 
(0, 1) be a sequence such that X^^Li 7 7« < °°- By Lemma ICTI there exists a null 
sequence of positive reals {k„} such that, for each fixed n, 

(6.3) liminf P q [ sup sup If (t) (d£ n ) < - ) > 1 - r\n- 

< i-° a \te[o,T]Cec nj 

For each n £ N, let M„ = {£ e M : sup CeC £(<9£") < 1/n}. If a sequence 
{£i} C M„ converges weakly to £, then for each open set <9£", the Portmanteau 
theorem yields 

£(0£») <lim sup UdT )<-■ 

i — >oo rt 

Thus, £ G M„ and M„ is closed. By definition of the Skorohod J\ -topology, the 
set T>1 = {((■) e D ([0,oo),M) : ({t) e M„ for all t G [0,T]} is also closed. Apply 
the Portmanteau theorem and l|6.3|l to obtain 



P (Z(-) G D£) > liminf P« (Z q {-) eDl)>l 



q — >oo 

By the Borel-Cantelli lemma, 



Vn 



\fc=l n=k J 

Thus, there exists a finite random variable N such that, almost surely, 

(6.4) sup sup Z(t)(d% n ) <-, for all n > N. 

te[o,T]Cec n 
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Since dc C d c n for all C G C and n G N, conclude that almost surely, 



sup sup Z(t) (dc) = 0. 

t£[Q,T}C£C 



□ 



We now establish property (JnJ. Recall that Z(t) — (l,Z(t)) for all t > 0, and 
S(u, v) = fj 1/Z(s) ds for all v > u > 0. 

Theorem 6.3. Almost surely, the process Z(-) satisfies 

(6.5) Z(t) (A) = Z(0) {A + {S(0, t), t)) + A f #(A+ (S(s, t),t- a)) ds, 

Jo 
/or allt>0 and A £ S3 . 

Proof. Let T > 0. It suffices to show that almost surely, (|6.5(l holds for all t £ [0, T] 
and all 4e J. For each r £1Z, define the random variable 

(6.6) X r T = sup sup £ r (M)(C) - \ r (t - s)^ r (C*) 

ceco<s<*<T 

pg 

By Lemma f5. II X^ — ► as q — > oo. Since the limit is deterministic, this conver- 
gence is joint with the convergence Z (•) =^ Z(-). Using the Skorohod represen- 
tation theorem, assume without loss of generality that {Z (-),Xj,} and Z(-) are 
defined on a common probability space such that 

(6.7) (Z q (-),X^) -> (Z(-),0), almost surely. 

The conclusions of Lemmas I6.ll and I6.2I hold almost surely as well. Assume for 
the remainder of the proof that all random objects are evaluated on the event of 
probability one such that Z(-) is continuous, and such that (|6.1|) . I|6.2|l and (|6.7I) 
hold. 

Fix t € [0, T] and C G C. An extension to all Borel sets A G SB will be made at 
the end. For each q, (|5.1U|) yields 

_____ 1 lE"{t) 

(6.8) z\t)(C)=Z q (Q)(c+(S q {0,t),t))+- _T l c (B q t (t),D q (t)] 

i—l 

We will obtain (|6.5() from l|6.8|l by letting q — > oo. The convergence in the first 
component of Ij6.7[l is in the Skorohod Ji-topology on D([0, oo),Mi). However, 
since Z(-) is continuous, 

(6.9) Z q {s) -^ Z{s), for all s G [0,t]. 
Since Z q (-) = (l,Z q (-)) and Z(-) = (1,Z(-)), this implies that 



= 0. 
t 



(6.10) lim Z H {-)-Z{-) 

q — >oo 

For all t > v > u > 0, (|6.1() implies that inf se r U]t ,i Z(s) > 0, and so the bounded 
convergence theorem yields 

— f v 1 
lim S (u, v) = lim / == ds 

9—°° < i^°°Ju Z (s) 

(6.1D _ r i ds 



u Z(s) 
S(u,v). 
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If Z(Q) 7^ 0, then Q6.11JI holds for u = as well, because then inf s6 [o^] Z(s) > 0. If 
Z(0) = 0, then 5(0, v) — oo and 5 (0, v) — > oo as g — > oo. 

Suppose that Z(0) ^ and let e > 0. By (|6.11|) . there exists ag e GK such that 
S 9 (0,t) e ((5(0, t) -e) + ,5(0,t) +e) for g > g e . Deduce from the shape of the set 
C, (JS2J, and (lOJ that 

limsu P Z>) (c+(S q (0,t),t)) < Z(0) (C+((S(0,t)-e) + ,t)) , 

(7 — UTO > ' 



lim inf Z 



By (fBTSf) . letting e -> yields 



(0) (C + (5 9 (0, i), t)) > Z(0) (C + (5(0, t) + e, *))) 



(6.12) lim Z q (0) (C + (5>, t),tj) = Z(0) (C + (5(0, t),t)) . 

If Z(0) = 0, then (|6.12l) holds trivially because the left side is bounded above by 
Um g ^ oo (l,Z 9 (0)) = by (|fTTU|l . Combining with (OJ and |JE2J for Z g (t), implies 
that, as q — > oo, 

Z g (i)(C) -Z 9 (0) (<7 + (5 9 (0,i),t)) - Z(t)(C) -Z(0) (C + (S(0,t),i)) . 

Let I q denote the second right hand term in l|6.8|l . Let 6 > and let r\ € (0, £). 
Since 5 (s, i) is decreasing in s and 5(-, t) is continuous on [rj, t] , (|6.11|) implies that 
5 (•,<) — > 5(-,t) uniformly on [r),t]. That is, there exists q$ € 1Z such that 



(6.13) sup S q (s 1 t)-S(s 1 t) 

selv-t] 



< 6, for all q > qg. 



Let D#(C) = {C e C : tf(d c ) ^ 0}. Note that A?(C) is countable because •&{■ x K+) 
and i?(R + x •) are probability measures. Since Z(u) > for all u G [r),t], the 
function S(s,t) is strictly decreasing in s on [r],i\. Thus, 

Do(S) = {s e [n, t] : C + (5(s, t) ± 2<J, t - s) € A?(C)} 

is also countable. For each integer N > 1, let r\ = £q < i^ < • • • < v§j = t be 

a partition of [77, t] such that t? £ D^(S) for all j = 1, . . . , N — 1, and such that 
maxj<jv— 1 (tj+i — if) — * as iV — > 00. Then 

<?£'(,,) AT-i g£ 5 (tf +1 ) 

iq = \ E ^ (s?«),E?(t)) + E J E l+ c (B q dt)Mt) 

Note that the first right hand term is bounded above by £ (0,ry)(IR + ). Suppose 
that tf < U^q^ 1 < t¥ + i, for some q > qs, some j < N — 1, and some z 6 
{qff{ri) + 1, . . . , g£ ? (£)}. Then by 1031 . 

(6.14) 5(if +1 ,t) - <5 < 5 9 ([/iV\i) < 5(tf ,t) + 5. 

By definition, 

(B?(t),35?(*)) = (B?-S'(U? q - 1 ,t) 1 DIq- 1 -(t-U?q- 1 : 
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So for q > qs, (|6.14J) and the inequalities lc(- — 5, •) < lj(-, •) < lc(- + S, ■) yield 
1+ (B?(t),I3«(t)) > 1 C (B? - (5(tf ,t) + 25), DU' 1 (t tf )) ; 
1+ (B«(t),25«(t)) < lc (B? - (5(tf +1 ,t) - 25), DU- 1 - (t- tf +1 )) . 

This yields, for q > qs, 

JV-l <lE\t» +1 ) 

iq ^ E - E J c (^ - (*(*f . *) + 2 <^ AV 1 - (* - if)) ; 
^<S 9 (o,v)K) 

JV-l 9^(tf+i) 

+ E ^ E !c ( B * - Wf+i > *) - 2«), AV 1 - (t - *f +1 )) . 
Rewrite as 

JV-l 

I 9 > E^f >**+i) (C+ (5(tf ,t) + 2*,*-tf)) ; 

<«■"> 

I* < Z*(0, »,)(<) + J] ^(*f >*&i) (C + (5(tf +1) t) - 2«J,t- *f +1 )) . 

By l|tj.6fl and l|6.15fl . q > qs implies that 

JV-l 

i q >J2 ( A9 (*f+i - *f )^ 9 ( c + ( 5 (*f »*) + 2 <^ - *f )) " x t) ; 

JV-l 

I* < X"r, + X« T +Y, (A 9 (tf +1 - if )& (C + (S(tf +1 ,t) -25,t- if +1 )) + X«) . 
By JE3, and since if £ A^S) for all j = 1, . . . , N - 1, 

JV-l 



liminf 1* > A £ ($_! - if )tf (C + (S(tf, t) + 25,t- tf)) ; 

limsupl 9 < At? + A £ (tf+i - *?)# (C + (S(tf +1 ,t) -2S,t- if +1 )) . 
For s £ [77, i] such that s £ D#(S) the bounded convergence theorem implies that 

JV-l 

w 1 i i ? oc E v ,*? +1 )(^ ( c + (^f . *) + 25 ' * - *f )) 

(6.17) 



= tf(C+{S{s,t) + 25,t-sj); 

JV-l 



^ E V >*f +1 )( 5 ^ ( C + CS^+i.*) -25,4 - *f +1 )) 
= 0(C+(S(«,t)-2tf,t-a)). 



JV 

J=0 
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Thus, the convergence in Q6.17JI holds for almost every s € [rj,t). Let N — > oo in 
(|6.16|l and conclude from (|6.17(l and the bounded convergence theorem that 

liminfF > A / tf(C + (S(s,t) +25,t- s)) ds; 

9^°° Jr, 

(6-18) JV t 

limsup/ 9 < A// + A / i? (C + (S(s, t) -26,t-s)) ds. 

q^oo J n 

Let 6 — ♦ in H6.18J1 . Since Z?,j(C) is countable, both integrands in H6.18|l converge 
almost everywhere on [77, t] to $ (C + (^(s, t), t — s)). Thus, 

liminf /«> A / #(C+(S(s,t),t-s)) ds; 

<?^°° Jr, 

limsup/«< Xr] + X I $(C+{S(s,t),t-s)) ds. 
Let f] — > to conclude that 

lim I q = A / i9 (C + (S(s, t), « - a)) ds. 

This proves igSJl for all i 6 [0, T] and C eC. To extend to all 4eJ, let C be 

the set of A e 58 for which l|6.5jl holds. Observe that C is a A-system: M + e C 

2 

because R + e C; if {A„} C C satisfies A n | A, then ^4 e C; if A\ C A 2 are 

elements of C, then ^2 \ A\ € C. Observe also that C is a 7r-system: if Ci, C2 € C, 

then Ci n C2 eC. Since CcC and the cr-algebra generated by C is equal to £§, it 

follows that C = S3 by the Dynkin 7rA-theorem (see for example 0). □ 
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